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ABSTRACT 

We present a detailed characterization of the X-ray spectral properties of 761 type 1 
AGN, selected from a cross-correlation of the SDSS DR5 quasar catalogue and the 
incremental version of the second XMM-Newton serendipitous X-ray source catalogue 
2XMMi-DR2. The X-ray spectrum of each source is fit with models based on a simple 
power law to which additional cold absorption and/or soft excess features are added 
if an F-test at 99% significance requires them. The distribution of best-fitting photon 
indices, T, is fit with a Gaussian with mean (r) = 1.99 ± 0.01 and dispersion cr^p) = 
0.30±0.01, however this does not provide a good representation of the distribution due 
to sources with very flat or steep T values. A highly significant trend for decreasing 
r values with increasing 2-10 keV luminosity, Lx , is seen but only a weak trend with 
redshift is found. Intrinsic cold absorption is detected in ~ 4% of the sample and soft 
excess emission is detected in ~ 8%. These values are lower limits due to the detectabil- 
ity being limited by the quality of the spectra and we suggest the intrinsic values may 
be as high as ~ 25% and ~ 80% respectively. The levels of rest-frame absorption are 
higher than expected for type 1 objects (Njj = 10 21 — 10 23 cm -2 ) and the fraction of 
absorbed sources and the Nu values are not seen to vary with Lx or z. The average 
blackbody temperature used to model the soft excesses is (kT) = 0.17 ± 0.09 keV. 
This temperature is found to correlate with Lx but not the blackbody luminosity or 
the black hole mass which do correlate with each other. A strong correlation is found 
between the luminosities in the blackbody and power law components suggesting that 
a similar fraction is re-processed from the blackbody to the power law component for 
the entire luminosity range of objects. A positive correlation between T and the X-ray 
derived Eddington ratio is found for the sources whose mass was determined using 
the H/3 line, but a negative correlation is found if the CIV line is used. No correlation 
is found if the Mgll line is used. No significant correlations are found between the 
blackbody temperature, luminosity or black hole mass with Eddington ratio, despite 
a link between the power law and blackbody production being indicated. The sample 
includes 552 confirmed RQQ and 75 confirmed RLQ. The RLQ are found to have 
higher Lx values than their RQQ counterparts, suggesting an additional X-ray com- 
ponent, perhaps related to a jet, is present in these sources. This component may also 
be the cause of the flatter T values seen in RLQ. 
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1 INTRODUCTION 

It is widely accepted that Active Galactic Nuclei (AGN) are 
accreting supermassive (M > 10 6 Mq) black holes (S MBH) 
residing at the centre of galaxies |Lvnden-Bellll 1969T ) . They 
have bolometric luminosities which can be as high as ~ 10 8 



erg s -1 (e.g. I Woo fc Urrvl [200^ 1 . with the X-ray emission 
(accounting for typically ~10% of the total luminous out- 
put) thought to originate from the innermost regions of the 
AGN. By studying their X-ray spectra we hope to further 
our understanding of the physical conditions surrounding 
the black hole. 
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The earliest observations of AGN revealed that their X- 
ray spectra were remarkably similar, and could be described 
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by a simple power law (|Mushotzkv et alj Il980h . This pri- 
mary X-ray emission is believed to be due to inverse Comp- 
ton scattering of low energy UV photons from the accre- 
tion disc. These photons are up-scattered to X-ray ener- 
gies by interaction with the corona, a cloud of relativis- 
tic electrons found above the surface of the accretion disc 
IIHaardt fc Maraschj|l993h . In higher quality X-ray spectra, 
additional features can be distinguished. Below 1 keV, excess 
emission with respect to an extrapolation of the power law is 
some times detected l|Arnaud et al.l 1 19851 ; I Turner fc Pounds! 
Il989h . whose origin remains unknown. The standard expla- 
nation for this 'soft excess' is that it is thermal emission 
from the inner accretion disc and hence is usually modelled 
as a blackbody component. The high energy tail of thermal 
emission from a standard thin accretion disc around a SMBH 
accreting at the Eddington rate is expe cted to have a black- 
body temperature of ~ 0.02 — 0.04 keV jShakura fc Sunvaevl 
1 19731 ). however the measured bl ackbody temperatures ar e 
typically higher than this (e.g. iGierlihski fc Done] l2004h . 
Lower energies are also affected by photoelectric absorption. 
In some sources absorption levels are detected above that 
which is expected due to the Galactic component. This may 
be due to cold or ionized (warm) gas intrinsic to the AGN, 
or due to its host galaxy. At higher energies (10 — 300 keV), 
a hump may also be detected due to the Compton reflec- 
tion of the incident power law continuum photons off the 
accretion disc. Despite being detected in the X-ray spectra 
of Seyferts, the evidence for this component in quasars is 
poor, since their emission is largely dominated by the lu- 
minous power law continuum. It has been suggested that 
the reflection component could also be intrinsically weaker 
in high luminosity quasars since the inner accretion disc is 
more likely to become ionized, thus reducing the amount 
of neutral matter rem aining, which is requir ed to create a 
reflection feature (e.g. iMushotzkv et aLlll993l ) 

The power law emission can be described by Pe oc E~ v , 
where the spectral (pho ton) index, F, was origi nally found 
to be approximately 1.7 (jTurner fc Poundsll 19891 ). As higher 
quality spectra became available the more complex features, 
including reflection , were discovered in addition to this un- 
derlying power law (I Pounds et al .1 19901 ). When the reflection 
component is modelled separately , the spectral index of th e 
power law is adjusted to V ~ 1.9 (|Nandra fc Poundsll 1994h . 
Many subsequent studies have attempted to constrain the 
intrinsic value of this photon index, and determine how var- 
ied this value is between different objects. 

In general, most X-r ay studies of AGN find an aver- 
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Dai et alj 12004 : iReeves fc Turner! 120001 ). However, different 
average values are reported for samples created by different 
satellites because of the selection effects occurring due to the 
different band passes. For example, samples selected in hard 
X-rays by INTEGRAL or Swift yield generally lower average 
spect ral indices of ~ 1.7 (|Molina et al.l [20091 : IWinter et al.l 
l2009h . likely because there is a bias towards detecting flat- 
ter sources at higher X-ray energies. The V value reported 
is dependent upon the energy range it is considered over, 



particularly if some spectral components are not properly 
modelled. The presence of a warm absorber or unmodelled 
reflection componen t can produce a slop e flatter than the 
true intrinsic value. iBianchi et all l|2009l ) suggest that the 
difference in photon index found between different popula- 
tions of AGN could be due to varying amounts of Compton 
reflection present in the sources. For their sample, quasars 
were found to have steeper power law slopes (F — 1.80±0.05) 
than Seyferts (F = 1.66 ± 0.05), which may be indicative of 
a larger reflection component present in the Seyfert pop- 
ulation. When a reflection component was included in the 
models, the F values of both populations increased, but more 
so for the Seyferts. 

Different T values in other populations of AGN are 
likely due to real physical differences. For example, ra- 
dio loud quasars (RLQ) are known to have flatter power 
law slopes than t h eir ra dio quiet counterparts (RQQ) (e.g. 
IReeves fc Turner! (|2000! ) RLQ: T = 1.66 ± 0.04, RQQ: 
r = 1.89 ± 0.05), likely due to their X-ray spectra being 
contaminated by synchrot ron emission from the ir power- 
ful relativistic radio jets (jSambruna et al] Il999l ). Narrow 
Line Seyfert Is (NLS1) are also observed to have steeper 
powe r law slopes than the ir broad- line counterparts (BLS1) 
(e.g. iBrandt et al] (|l997i ) NLS1: F = 2.15 ± 0.05, BLS1: 
r = 1.87±0.04), which could be due to these objects accret- 
ing at high rates. For higher mass accretion rates, a stronger 
soft X-ray/UV component from the disc is created, which 
in turn causes more Compton cooling, producing a steeper 
X-ray spectrum. This scenario was offered as an explana- 
tion for the NLS1 R E 1034+39, which h as an unusually soft 
power law of r ~ 2.6 (|Pounds et al.ll 19951 ). The spectral state 
observed was noted to be similar to the high/soft state ob- 
served in galactic black hole binaries that are accreting close 
to the Eddington limit (|Remillard fc McClintockll2006l ). 

Despite being widely investigated, it is still not clear 
whether there is any intrinsic variation of F with X-ray 
luminosity or redshift. It is natural to assume that the 
X-ray production is related to the black hole mass, ac- 
cretion rate and l uminosity, howev er this dependence is 
poorly understood. iDai et al.l l|2004l ) first reported a cor- 
relation with luminosity that was positive, implying that 
the X-ray spectral slopes were steeper in higher luminos- 
ity sou rces, however later studies eith er rep ort no corre- 
lation (IWinter et al.l 120091; ITozzi et al.l 120061; IMateos et~aH 
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ray spectral slope flattens as the X-ray luminosity of the 
sources increases. If a trend with redshift was found it 
may indicate that the accretion mechanism was different 
for objects in the past. Many studies of larg e samples of 
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which agrees with some results of studies specifi- 



cally targeting a few high redshift objects and measuring 
F values consistent with those for lower redshift objects 
(|Vignali et alj|2005l ; Ishemmer et aT1l2005l ). However, some 
report F flattening for higher redshift sources (|Kellv et al.l 
l2007l : iBechtold et alJliooi ; IVignali et al.lll999l) 

The Unified Model of AGN (|Antonuccil 19931 ) states that 
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different classes of AGN are examples of the same physi- 
cal object but viewed at different angles. In type 1 objects 
we are able to observe the inner most region, whereas in 
type 2 objects this is shielded from view by a dusty torus. 
Therefore type 1 AGN are expected to be unabsorbed in 
X-rays and hence require no intrinsic absorption compo- 
nents in the modelling of their spectra. However, previous 
studies have found bo th type 2 objects which appear to 
be X-ray unabsorbed (|Panessa fc Bassanil 120021 ') and sam- 
ples of type 1 AGN with significant fractions of absorbed 
sources (|Mateos et al.ll2010l ; lYoung et al .120091 ; IMateos" et al.l 
l2005bl l. These results are not reconciled by the standard ori- 
entation based Unified Model. 

Soft excess components are found to be common fea- 
tures in the X-ray spectra of type 1 AGN. Originally 
thought to be present in 30% of hard X-ray selected Seyferts 
|Turner fc Pounds! [l989l ), more recent estimates give ~ 7% 
(e.g. IMateos et al]|2010l l. However the origin of this compo- 
nent is still a matter of debate. It was originally interpreted 
as the hard tail of the 'Big Blue Bump' seen in the UV due 
to the accretion disc, but is now thought to be either an arte- 
fact of ionized absorption JCierliriski fc Pond 20041') or io n- 
ized reflection (|Ross fc Fabianll2005l ; ICrummv et alj|2006h . 

One of the key questions in AGN research is to explain 
why different classes of AGN have different spectral prop- 
erties. The best way to do this is though a stat istical study 
of a large number of AGN. iBianchi et all (|2009l ') investigate 
the X-ray spectral properties of ~ 160 AGN from targeted 
XMM-Newton observations. This sample has the advantage 
of high signal-to-noise (S/N) X-ray spectra, however it is also 
biased towards sources at low redshift (90% of their sources 
are at z < 1). Since XMM-Newton targets have been used, 
the sample contains objects observed for a number of differ- 
ent reaso ns, and hence does n ot represent a uniform sample 
of AGN. IMateos et alj (|2010T ) presented an X-ray spectral 
analysis of ~ 500 type 1 AGN from the XMM-Newton Wide 
Angle Survey (XWAS). The overall broad-band properties 
and their dependences on X-ray luminosity and redshift were 
investigated, offering a complete picture of the X-ray spec- 
tral properti e s of a n X-ray selected sample of type 1 AGN. 
lYoung et all (|2009l ) present an optically selected sample of 
~ 500 objects by cross-co rrelating the Sloan Digital Sky 
Survey (SPSS; lYorklbOOOT ) and the XMM-Newton archive. 
They fit only simple models to their spectra and include only 
~ 300 of their sources in their investigation of T with Lx ■ It 
is possible that by restricting their sample to the higher S/N 
sources, they have introduced a bias towards lower redshift 
objects, or the most X-ray luminous higher redshift sources. 

In this study we analyze the X-ray spectra of the largest 
sample of type 1 AGN to date (~ 750). They are optically 
selected from SDSS and also have observations taken by 
XMM-Newton. We fit a range of spectral models to each of 
our sources, over an energy range of 0.5-12.0 keV, including 
ones to model any soft excess components in order to gain 
insight into their nature. We choose the best-fitting model 
for each source and constrain the average power law slope 
and intrinsic dispersion. We also investigate any absorption 
found in the sources. Our sample covers a large range in 
redshift, with sources up to z ~ 5, and covers around 4 or- 
ders of magnitude in Lx- This allows us to investigate the 
dependencies of Y with X-ray luminosity, Lx, and redshift, 
z. 



This paper is organized as follows. In Section [2] we de- 
scribe the properties of the two initial catalogues, DR5QSO 
and 2XMMi and the extraction of the X-ray spectra and the 
subsequent spectral fitting. Sections 13.11 13.21 and 13.31 con- 
sider the overall distribution of power law slope values and 
the correlations of T with redshift and luminosity. Section 
3.41 explores the radio properties of the sample. Section 13.51 
considers the properties of the sources which required an 
absorbed power law model, whilst those requiring a soft ex- 
cess component are discussed in Section [3.61 In Section 13.71 
we use black hole mass estimates from the literature to in- 
vestigate how the accretion properties affect the observed 
X-ray emission. A discussion of our results and conclusions 
are presented in Sections [4] and [5] respectively. 

In this paper we assume a cosmology from the WMAP 
resu lts where Hp=70 km s _1 Mpc -1 , Q, M = 0.3 and = 
0.7 jSpergel et alj|2003h . 



2 PROCEDURE 
2.1 Data Sample 

The sample of objects for this study was created by a posi- 
tional cross-correlation of the SPSS DR5 quasar ca t alogu e 
(|Schneider et all 120071 ') and 2XMMi (| Watson et all l2009h . 
using a 10" matching radius. This initial sample consisted 
of 2013 matches, corresponding to 1281 unique objects. 

The PR5 quasar catalogue is produced fr om the fifth 
data release of SDSS (| Adelman-McCarthyfeOOTT ) . which pro- 
vides imaging in 5 broad optical bands (ugriz). Quasar can- 
didates are primarily selected for spectroscopic follow up 
based on their location in multidimensional SDSS colour 
space, however others are selected if they lie within 2.0" 
of a FIRST (Faint Images of the Radio Sky at Twenty- 
Centimeters; [Becke^et^ai][l99^) radio source, introducing a 
possible radio bias into the catalogue. 6.6% of our final sam- 
ple of sources were selected for follow up due to a close match 
with a FIRST detection, however all but 2 of these sources 
would have been included automatically by the colour selec- 
tion. The catalogue includes objects with an absolute i band 
magnitude brighter than —22, an apparent i band PSF mag- 
nitude fainter than 15, a reliable redshift and at least one 
emission line in their spectra with a FWHM > 1000 fans" . 
77,429 quasars are included with redshifts from 0.08 to 5.41, 
covering a sk y area of about 5740 deg 2 . 

2XMMi jWatson et alj|2009f ) is the incremental second 
serendipitous X-ray source catalogue, consisting of both the 
targets and non-target source detections that are recorded in 
the large (30') field of view of XMM-Newton' s EPIC (Euro- 
pean Photon Imaging Camera; lTurner|[200ll ; IStriider et al.l 
l200ll ) X-ray cameras. The catalogue contains information 
from 4117 pointed observations covering a non-overlapping 
sky area of ~ 420 deg 2 . It contains ~ 290,000 detections, 
which due to field duplication by multiple observations, cor- 
responds to ~ 220, 000 unique sources. Since observations 
with XMM-Newton can be requested for a variety of rea- 
sons, including the target detections in our sample could 
add bias towards 'interesting' objects. Our final sample in- 
cludes 62 sources which were the target of an observation 
(8% of the sample), which are naturally the ones with the 
highest number of counts and the sources fit with the most 
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complex models. Since our analysis deals primarily with the 
distribution of Y values and trends of Y with various quan- 
tities, we investigate whether the target sources have a sig- 
nificantly different Y distribution compared with the sample 
of serendipitously detected sources. A Kolmogorov-Smirnov 
(KS) test finds them to be consistent with each other and 
therefore the use of the target detections in our data sample 
should not affect our overall findings. 

We determine the radio loudness of the objects in our 
sample using the definition Rl = Fr/Fo where Fr is the 
flux at 5 GHz and Fo is the flux at 4400 A. Any sources for 
which Rl > 10 are taken to be radio loud (|Kellermann et all 
I1989D . To do this we cross-correl ated our X-ray sample with 
the FIRST (|Becker et al.lll995h catalogue which covers a 
similar area to SDSS at a frequency of 1.4GHz to a sen- 
sitivity of ~ 1 mjy. 128 quasars had a match with a FIRST 
source, 44 fell outside the FIRST area, and for the remain- 
ing sources 5 x RMS value at the source position was used 
as an upper limit. The flux measured by FIRST was extrap- 
olated to a flux at 5 GHz by assuming the radio emission 
is described by a power law with spectral (energy) index 
a = —0.8. The flux at 4400 A was determined from a linear 
interpolation of the optical fluxes in each of the SDSS ugriz 
bands. 552 (72%) were found to be radio quiet quasars, 75 
(10%) radio loud quasars, 90 (12%) were sources whose up- 
per limits are too high to determine their radio loudness 
with confidence and 44 (6%) sources lay outside the FIRST 
coverage area and therefore a radio loudness could not be 
determined. If all the sources with radio upper limits are as- 
sumed to be either all radio quiet or all radio loud, we find 
that the radio loud fraction lies between ~10% and ~22%, 
respectivel y, in broad agreement w ith previously reported 
values (e.g. iKellermann et al.|[l989l ). 

2.2 Extraction of Spectral Products 

The X-ray spectra of all 2013 detections were extracted using 
an automatic pipeline incorporating tasks from the Science 
Analysis Software v6.6.0 (SAS3), specifically designed for 
the reduction of XMM-Newton data. In particular the task 
eregionanalyse was used to select the extraction regions 
so as to optimize the S/N. Circular regions with an average 
radius of ~ 17"were used for the source spectra. A circular 
background region with a radius of 50" was automatically 
chosen in the same CCD and close to the source, but ex- 
cluding any other nearby sources. Auxiliary response files 
and redistribution matrices were generated using arfgen 
and rmf gen, respectively. 

Each detection was subjected to a manual screening 
process in which any that were severely affected by Out of 
Time events or contamination from a nearby bright source 
were removed from the sample. Sources falling into chip gaps 
or partially off the edge of the field of view were also re- 
moved. Multiple detections of the same source were merged 
into a single spectrum, with the data from the two MOS 
cameras being combined, but a separate MOS and pn spec- 
trum for each source generated, assuming that data from 
both detectors was available. For 278 sources data from only 



one detector was available or hadn't been rejected in the 
screening process. For these sources, spectral fitting was car- 
ried out on the single spectrum available. 6% of the original 
sample of unique objects were removed in the screening pro- 
cess, leaving a sample of 1201 AGN. Objects with < 75 total 
counts were excluded, leaving 761 sources. 

2.3 Spectral Fitting 

The spectra were binned using the grppha Ftoo0, to ensure 
that Chi squared (\ 2 ) statistics could be used. The min- 
imum number of counts per bin was varied depending on 
the S/N of the spectra, with 15 being used for sources with 
~ 100 counts, increasing up to 30 counts per bin for sources 
with a few thousand counts. The spectra were then fit with 
phys ically motivated models using XSPEC vll.3.2 |Arnaudl 
Il99rj ) over the energy range 0.5-12.0 keV. Where available, 
the MOS and pn data were fit simultaneously with the same 
parameters and a freely varying constant was added to the 
model to account for calibr ation offsets between the two 
cameras |Mateos et al.l [20091 ) . Each model included an ab- 
sorption component fixed at the value of the Galactic column 
density in the direction of the source, as determined from 
the nH FtooP^ which uses the HI map of lDickev fc Lockmanl 
|l99d ). 

The 761 sources with more than 75 counts were fit 
with both a simple power law ('po') and an absorbed power 
law ('apo'), in which Y and an intrinsic Nh component at 
the redshift of the source were allowed to vary freely. For 
680 sources with more than 100 counts, models including 
a blackbody component were also considered. This com- 
ponent was introduced to model any soft excess, and was 
added to both the simple power law and absorbed power 
law models. A summary of the different models considered, 
along with their XSPEC terminology can be found in Ta- 
ble [1] Absor ption was m o delled using phabs with abun- 
dances from IWilms et al.l (120001) and c ross-sections from 
iBalucinska-Church fc McCammonl (| 19921 ). 

The best-fitting model for each source was chosen via 
an F-test at 99% significance. The best fit was assumed to be 
a simple power law unless an F-test determined that an ad- 
ditional intrinsic absorption component or blackbody com- 
ponent was required. These sources now assumed to be best 
fit with the 'apo' or 'po+bb' models were then tested to see 
if the best fitting model was in fact 'apo+bb' i.e. both com- 
ponents were required. The number of sources best fit with 
each model can also be found in Table [T] Example spectra of 
sources best-fit with different models are shown in Figure [T] 

While the F-test determined the best fit to the data, a 
model was deemed to be a good fit to the data if the null 
hypothesis probability (HO) was > 1%. The majority (96%) 
of the sample were well fit by the model chosen by the F- 
test and the numbers of good fits for each type of model 
can be seen in Table [1] For the sources with high numbers 
of counts (i.e. > 1000 degrees of freedom) the models we 
fit are not sufficient to deal with the complex ities shown to 
exist from studies of bright, nearby AGN (e.g. IPounds et al.l 
l2004h . and therefore HO < 1%. For a good fit we expect 
X 2 /v ~ 1, which is shown in Figure [5] and compared to the 



1 The description and documentation arc available online at 
http: //xmm2.csac.csa.mt / sas/ 



2 http:/ /heasarc. gsfc.nasa.gov/docs/softwarc/ftools/ 
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Figure 1. This figure includes the raw (left) and unfolded (right) spectra for 3 of the sources in the sample which are each best fit 
with different models. 2XMM J092633. 4+304535 (top row) is best fit with the simple power law model (po), where T = 1.98 ± 0.07. 
2XMM J085724. 0+090349 (middle row) is best fit with the absorbed power law model (apo), where F = 1.85 ± 0.36 and N H 
(9.8 + 3.6) x 10 22 cm~ s 
T = 1.90 + 0.09 and kT 

left panels, the MOS data is shown in black and the pn in red. In the bottom right panel, the dashed line indicates the contribution of 
the blackbody component to the overall spectral model. 



2 . 2XMM J154316. 4+540526 (bottom row) is best fit with a power law model and a soft excess (po+bb), where 
0.10 + 0.02 keV. The sources include (from top to bottom) 630, 129 and 1131 total (MOS+pn) counts. In the 
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Table 1. This describes the different models used to fit the X-ray spectra of the sources. The spectral components included in each 
particular model arc listed in column 2. Galactic absorption is also included in each, modelled by the additional phabs component. 
The 4th column lists the total number of sources in the sample which were best fit by that particular model as chosen by an F-test 
at 99% significance. The final column lists the number of those sources which are also well fit by the best-fitting model i.e. HO > 1%. 
We note that the blackbody component is only considered in spectra with > 100 counts (see Section 13.611 



Model 


Components 


XSPEC terminology 


No. of sources 


No. of good fits 
HO > 1% 


po 


Power Law 


phabs*po 


672 


650 


(89%) 


apo 


Power Law + intrinsic absorption 


phabs*zpha*po 


2 gO 


28 


(4%) 


po+bb 


Power Law + soft excess 


phabs* (zbbody+po) 


55* 


52 


(7%) 


apo+bb 


Power Law + intrinsic absorption - 


- soft excess phabs* (zbbody+zpha*po) 


t 


■1 


«1%) 








761 


734 





1 source shows a Fe Kct line and is fit with the model pha*(zpha*po+zgauss) 

2 sources show Fe Kct lines and are fit with the model pha*(zgauss+zbbody-|-zpha*po) 



actual quality of the fits. The region in which the best fit 
model is considered correct at the 99% confidence level is 
shown between the red dashed lines. 

The spectral fits of each source, particularly those with 
a poor original fit, were visually inspected and revealed 3 
sources with evidence for the Fe Ka emission line. A Gaus- 
sian component, zgauss, was manually added to the model 
for these sources. We note that evidence for this component 
likely exists in more of the source spectra, although a sys- 
tematic search for this component is beyond the scope of this 
paper. The 3 Gaussian components have energies (E) and 
equivalent widths (a) as follows: E = 6.72 ± 0.002 keV and 
a = 0.07±0.01 keV, E = 6.67±0.03 keV and a = 0.22±0.03 
keV, E = 6.40 ± 0.05 keV and a = 0.13 ± 0.02 keV, which 
are typical of those found for nearby, bright AGN. 

Fluxes and rest frame luminosities for the sources are 
computed from the best fit model which were corrected for 
Galactic absorption, and in the case of the luminosities, in- 
trinsic absorption where necessary. The final sample covers 
a large range in redshift (0.11 to 5.41) and 2-10 keV X-ray 
luminosity (10 43 — 10 46 erg s~ 4 ), better enabling us to inves- 
tigate the dependencies of T on these quantities. The redshift 
and counts distribution can be seen in Figure [3] (top) and 
the luminosity-redshift dependence can be seen in Figure [3] 
(bottom). The properties of our sample are consistent with 
those of the entire DR5 quasar catalogue in terms of both 
the average redshift and average absolute i band magnitude. 
More details will be presented in Hutton et al. (in prep). 

As expected, the majority of our sources, ~ 90%, are 
best fit with the simple power law model and do not require 
any additional spectral components. Figure|3]shows the frac- 
tion of sources that are detected with either an absorption 
component or a soft excess component as a function of the 
number of counts in the spectra. At the highest count lev- 
els approximately 25% of the sources require an absorption 
component while over 80% require a soft excess. Both frac- 
tions fall with decreasing counts. As this could be due to 
our reducing sensitivity to these extra components we have 
performed simulations of the detectability of typical values 
of these as a function of counts (Scott et al. in prep). We 
find that these simulations reproduce almost exactly the be- 
haviour observed in Figure [4j suggesting that the level of 



000 



100 - 



x 




10 - 



00 



000 



V 



Figure 2. Plot of Chi squared value (x 2 ) vs Degrees of freedom 
(is) for each of the X-ray spectra. The solid line indicates the 
location of x 2 / u = !• The dashed red lines represent the locus 
outside which we would expect to observe \ 2 values only 1% of 
the time if the fit is correct. 



required soft excess and absorption within the overall type 
1 AGN population could be as high as the 80% and 25% 
found in the highest count range. 

We have also conducted a joint spectral fitting of 
sources in narrow redshift bins to see if the additional spec- 
tral features are detected in a combination of lower count 
(< 500) spectra which did not previously require them in 
the fit (i.e. were best-fit with the simple power law model). 
We did this by fixing F = 2.0 for each source and fit the 
apo (or po+bb) model to the spectra, to find a single Nh 
(or kT) value for all objects. We expect from Figure [4] that 
fitting spectra with ~ 7000 combined counts, will reveal the 
soft excess feature, as it should be detected 25 — 40% of the 
time. We find that for the lower redshift bins (z < 1), the 
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Figure 3. Plots showing the general properties of the final sam- 
ple of 761 objects with X-ray spectral fits. Top - The redshift 
and counts distribution. Bottom - The luminosity (2-10 keV) and 
redshift range of the sample. 



addition of a soft excess component, does result in a bet- 
ter fit to the combined spectra and the component has a 
temperature consistent with those reported in Section 13.61 
At ~ 7000 counts, we only expect 3 — 10% sources to be 
detected with an absorption component. Therefore in our 
sub-bins of ~ 30 sources, we do not expect to recover a sig- 
nificant absorption component. This is what we find, with 
no significant improvement of the fit upon the addition of 
an intrinsic Nh component. Further details of this analysis 
will be presented in Scott et al. (in prep). 



3 RESULTS 

3.1 Distribution of Power Law Slopes 

The distribution of best-fitting power law slope values, V, 
can be seen in Figure [5] It shows that most sources have 
a power law slope value ~ 2, with the arithmetic average 
value being V = 1.97 ± 0.01. The typical (68%) error on a 
r measurement is AF = 0.13, significantly smaller than the 
spread of F values, indicating that the observed dispersion 
of values is truly intrinsic. The panel showing the enlarged 
distribution indicates that there are some sources with ex- 
tremely flat (F < 1.4) or steep (T > 2.6) power law slopes. 

Such extreme values may be due to sources having spec- 
tra with low numbers of counts and hence the quality of the 
fit is poor. Figure [5] shows the dependence of the derived V 
value on the quality of the spectra, and it can be seen that 
many of the extremely flat or steep sources are indeed those 



CD 

cn 
o 

c 

CD 

(J 

CD 

Cl 

TJ 
CD 

O 
CD 

CD 
Q 



1 



Absorption 

Soft Excess 

1% Spurious Level 



10- 5 



Counts 



Figure 4. This plot shows how the percentage of sources that 
are detected with an absorption or soft excess component varies 
depending on the number of counts that are available in the spec- 
tra. 



with lower quality spectra. This is likely due to the presence 
of intrinsic absorption or a soft excess component which are 
not detected in the spectra with enough significance. How- 
ever, we note that there are sources with high counts for 
which we do not believe this to be the case and therefore 
do have intrinsically flat or steep V values. A more detailed 
discussion of these sources is beyond the scope of this paper. 

We assume that we can represent the distribution of V 
values with a Gaussian function of mean (F) and dispersion 
<T(r) and obtain simultaneous estimates of these values using 
a x 2 minimization technique. Figure [5] (bottom) shows the 
la, 2a and 3a confidence contours, along with the derived 
best-fitting values (r) = 1.99 ±0.01 and a {r) = 0.30 ±0.01. 
The best-fitting Gaussian can be seen overplotted on Fig- 
ure [5] (top). 

This analysis does not take into account the individual 
errors on each of the F estimates, so here we investigate what 
contribution these make to the total dispersion we measure. 
This is done by assuming that the observed dispersion can be 
expressed as the sum of the intrinsic dispersion and the indi- 
vidual errors on F added in quadrature, i.e. a 2 = a 2 nt + a 2 rr . 
We then calculate \ 2 f° r a range of different values of ai nt , 
assuming the expected F value for each point to be equal to 
the average value for the entire sample and the error on each 
point to be the 68% error derived from its spectral fit. The 
results of these fits can be seen in Figure [7] As expected, 
we find a 'poor' fit for low intrinsic dispersion, while the 
fit becomes 'too good' for much higher levels. The point at 
which reduced x 2 becomes equal to 1 is taken as the best 
fitting value and we use the null hypothesis probabilities to 
determine the 68% confidence interval on this value. We find 
aint = 0.31 ± 0.01, which is the same result that we obtain 
from our single Gaussian fit above. We can therefore say 
that the individual errors on T contribute very little to the 
overall dispersion. This is not entirely unexpected since the 
typical 68% error on T is Ar ~ 0.13, which will contribute 
much less than an intrinsic dispersion of 0.30 when added in 
quadrature. 

Our values of F agree with those found by previous 
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Figure 6. This shows the measured T values against the total 
number of counts in each of the X-ray spectra for which HO > 1%. 
The error bars have been omitted from this plot for reasons of 
clarity but an error of typical size can be seen in the top right 
hand corner. 
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Figure 7. This shows how the reduced \ 2 of the fit varies when 
increasing levels of intrinsic dispersion are added in quadrature to 
the individual errors on the T values for each point in the sample. 



Figure 5. The histogram (top) shows the distribution of best- 
fitting photon indices where the required model gives an accept- 
able fit to the data. The smooth line is the best-fitting Gaussian 
distribution, as determined by x 2 minimization and being de- 
scribed by (r) and 0"(r> ■ Also included is a plot of the residuals, 
comparing the data to the model, and an enlarged version of the 
lower regions of the histogram to better highlight the presence 
of the extreme sources. The bottom figure shows the best-fitting 
values, indicated by the cross, and the la, 2a and 3a confidence 
contours for the derived values of (T) and a^ r y 



studies, includi ng the XWAS sam ple, T = 1.96 ± 0.02, 
a = 0.27+Q 02 (|Mateos et alj|201fj ) a nd the XMM-Newto n 
COSMOS sample T = 2.06 ± 0.08 (jMainieri et alJ l2007t ). 
lYoung et all (|2009T ). using a subset of objects from the same 
data presented here, find T — 1.91 ±0.08 and that their typ- 
ical 68% errors, AF = 0.15 result in an intrinsic dispersion 
of a = 0.37. 

When fitting a single Gaussian to the T distribution, 
a null hypothesis probability value of p=0.43% was found. 
Since this is < 1%, it is assumed that the model does 
not provide an acceptable fit to the data. A two Gaussian 
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model gives a better fit to the data (null hypothesis prob- 
ability, p=15.4%) and an F-test comparing the two Gaus- 
sian and one Gaussian models suggests that the inclusion of 
the second Gaussian is statistically valid (F-test probabil- 
ity of 99%). In this fit the average value of the two Gaus- 
sians, which were kept equal, was p = 1.98 ± 0.01, and the 
best-fitting dispersions of the Gaussians were found to be 
o\ = 0.27 ± 0.01 and a 2 = 0.59 ± 0.04. The fit can be seen 
in Figure [8] We do not suggest that the 2 Gaussian model 
should be taken as a physically motivated fit to the Y distri- 
bution since it models both the fiat and steep sources as a 
single 'extreme' population, when in reality they are likely 
to be different sub-samples. However since a multiple Gaus- 
sian model can produce an acceptable fit, where the single 
Gaussian is inadequate, it does highlight the existence of 
extreme sources in the distribution. We test whether these 
are simply those with the lowest numbers of counts and find 
that a single Gaussian will provide an acceptable fit to the 
r distribution only when sources with less than 300 counts 
are excluded, at which point the sample size is reduced to 
385. At this count limit, we find that a = 3.7 x AF, where 
AF has fallen to 0.07, indicating that the dispersion is still 
significantly larger than the typical error. 

3.2 Redshift Dependence of the Power Law Slope 

In this Section we consider whether our data shows any evi- 
dence for an intrinsic variation of the spectral index, F, as a 
function of redshift. This possibility has been w idely investi- 
gated, but no clear co nsensus has been reac hed dYoung et al.l 
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eludes the 734 sources for which the best-fitting model pro- 
vides an acceptable fit to the data i.e. HO > 1%. The sources 
are plotted with different colours/symbols according to their 
model. 

We calculate Spearman's rank correlation coefficient, 
p, for both the entire population of 734 sources, and the 
subset of 650 sources whose best fit is a simple power law 
model (po). Considering the entire population we find a ~ 
3cr anticorrelation of V with z (p = —0.11, prob = 0.002). 
Considering only the sources fit by po, we find a tighter ~ 5a 
anticorrelation (p = -0.20, prob = 4 x 10 -7 ). 

Figure [10] shows binned versions of Figure [9] The top 
plot includes sources fit with all models, whilst the bottom 
plot only considers sources best fit with the simple power 
law model. This is considered because we find that a large 
percentage of sources in the redshift bins z < 0.6 are those 
fit with the po+bb model, whereas the bins at higher red- 
shift are mostly sources that are fit with the simple power 
law model. Whilst the inclusion of the soft excess compo- 
nent should account for the spectral complexity present and 
give r values that are consistent with those seen in the sim- 
ple power law sources, this is not found for our sample. The 
average F value for the sources fit with the po+bb model is 
significantly lower than that of the entire sample as deter- 
mined by a KS test. 

Narrow redshift bins (d£=0.2) are used up to z — 2.6. 
Above this point, we bin so as to ensure at least 8 sources 











30 




20 


(n 




D 


10 






TJ 





'</) 




<D - 


10 


a: 




20 




30 


cn 



30 


o 




O 


20 


in 




o 


10 


6 





z 





2 

r 



Figure 8. The red, dashed and blue, dotted lines show the two 
constituent Gaussians, while the black, solid line shows the overall 
model, produced by the addition of the two. 



per bin due to decreased numbers of sources, (or 6 in the 
case of the po only plot). A weighted mean value of F is 
calculated for each bin according to the equation: 



<r)=^P i xr i where P, = 



(1) 



which takes into account the individual errors (<7i) on 
the r estimates (r^). The errors on (r) are calculated as 
standard errors on the mean, a = a/yN, where: 



N 



(2) 



We fit the data with 2 models; a linear trendline where 
a gradient inconsistent with zero would indicate F values 
evolving with z, and a model where (F) was fixed at the value 
of the weighted mean for the sample to represent a non- 
evolving F. The models and corresponding null hypothesis 
probabilites (p- values) can be found in Table [2] 

For the sources fit with all types of spectral model, a 
non-evolving V gives a good fit to the data and the linear 
trend gives a gradient consistent with zero. This suggests no 
significant trend for decreasing V with increasing z. However, 
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Table 2. Listed are the best-fitting linear correlations between V and a or Lx for the binned data, found by x 2 minimization. The 
errors on the values are la. The null hypothesis probabilites (p-values) for both the linear models, and constant (r) models are 
listed in column 4, where p> 1% is taken to be an acceptable fit. 



Correlated Properties 
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Figure 9. This plot includes the 734 sources with H0> 1%. The 
different colours/symbols correspond to the best fitting model; 
red circles are po, green squares are apo, blue diamonds are 
po±bb and black triangles are apo±bb. The error bars shown 
at the top of the plot correspond to the average of each of the 1<t 
errors on the individual sources, calculated in dz=0.5 bins. 



when only the sources fit with a simple power law model are 
considered, the non-evolving V model does not provide a 
good fit to the data and the linear trend gives a gradient 
inconsistent with zero. This suggests that for this sub-set of 
sources, a significant, but slight, trend of F with z is present. 

We note here that our estimate of (r) for the entire sam- 
ple is lower than the arithmetic value determined without 
considering the individual errors on Y measurements. This 
does not seem to be due to any particular source with vastly 
underestimated errors, but a general effect created because 
the slopes of flatter sources are easier to constrain and there- 
fore have lower associated errors than the steeper sources. If 
we remove the top 2.5% of sources with weights above 0.01, 
we find the weighted mean value steepens to 1.93. 

The lowest redshift bin of Figure [10] has a large error in 
(r) due to a low number of sources with a large scatter of V 
values, which also causes the (F) value to change consider- 
ably between the 2 plots. Similarly the higher redshift bins 
also include low numbers of sources and are likely unrepre- 
sentative of the entire sample as the percentage of sources 
that were targets of the XMM-Newton observation is higher 



(for z > 2.6 : 25 ± 11%, for z < 2.6 : 7 ± 1%). Therefore 
we investigate whether the lowest bin and top bins strongly 
influence the results and find that the overall conclusions 
remain the same, i.e. there is no dependence of F on z when 
all sources are considered, but a marginal trend is apparent 
when just those sources best fit with a simple power law are 
considered. 

Since it is suggested that RLQ tend to have flatter spec- 
tra than RQQ, it is possible that an increased fraction of RL 
sources in high redshift bins may be a contributing factor to 
the flatter F values observed. However, all the sources in our 
sample above z — 2.6 are radio quiet using the definition of 
Rl from Section [2] 

We note that when lower count sources are removed, 
the results for the binned data remain similar. 

As our sample is flux limited, we are biased towards 
higher luminosity objects in higher redshift bins. We there- 
fore consider in Section \3 .3 1 whether we are observing flatter 
power law slopes in higher luminosity objects, rather than 
flatter power law slopes in higher redshift objects. 

As we find a trend of F with z, we investigate whether 
this could be due to the presence of a reflection component 
which may be detected in higher redshift sources where the 
EPIC bandpass probes up to higher rest frame energies. If 
this is not modelled separately, and just a simple power law 
model is used for the entire spectrum, its presence will make 
the slope of the power law appear flatter. This could explain 
the slight decrease in the power law slopes seen for sources at 
higher redshift. To test this hypothesis we fit simple power 
law models to spectra over the rest frame energy range 2.0- 
6.0 keV, with the upper limit chosen so as to exclude any 
contribution from such a reflection component and give a 
cleaner representation of the power law. The lower limit of 
2 keV is chosen to exclude any absorption component. 

Figure [TT] plots the power law slope value measured 
over the reduced range, 2.0-6.0 keV in the rest frame (r part ) 
against the slope measured over the full 0.5-12.0 keV ob- 
served frame (Ff u n), for a reduced sample of 128 sources 
previously best-fit with a simple power law, and including 
> 250 counts over the reduced range in each camera that 
is in use. They are divided into 3 broad redshift bins and 
the typical 68% error size is indicated to the left of the fig- 
ure. If we are removing an unmodelled reflection component 
when considering the spectral fits over the reduced range, 
we would expect r pa rt to be steeper than Ff u n and hence 
the sources should lie in the top left corner, above the r par t 
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Figure 10. These plots show the dependence of weighted T values 
on redshift. The sources are binned in redshift (dz=0.2), up until 
z = 2.6 and then after this, to ensure at least 8 sources per bin 
(or 6 for the po only plot). The y error bars are standard errors 
on the mean, a = <r/\/~N, where the standard deviation is based 
on deviations from the weighted mean. The x errors indicate the 
width of the bin. The solid horizontal lines indicate the weighted 
mean values of each sample and the dotted lines shows the linear 
fits to the data. The top plot includes data from 734 sources with 
good fits and a linear fit of (r) = (-0.04 ± 0.02)2 + (1.94 ± 0.03). 
The bottom plot includes data from 650 sources with good fits 
using a simple power law model and a linear fit of (T) = (— 0.07± 
0.02)2 + (1.99 ±0.03). 



= Tfuii line. We calculate the x statistic between the points 
and the null hypothesis line of r part = Ff u n. We find all 3 
redshift bins to be consistent with this line indicating that 
there is no significant difference between Tfuii and r part . This 
suggests no strong reflection component is present in these 
sources and therefore cannot be used to explain the flatten- 
ing of r with increasing z. 



3.3 Luminosity Dependence of the Power Law 
Slope 

In this Section we consider whether our data shows any evi- 
dence for an intrinsic variation of the spectral index, F, as a 
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Figure 11. This figure plots the power law slope value as de- 
termined from a spectral fit over the entire 0.5-12.0 keV detector 
energy range against the value determined from a spectral fit over 
the reduced 2.0-6.0 keV rest frame energy range for 128 sources 
as detailed in the text. They are divided into 3 broad redshift 
bins, red circles: z = 0.1 — 0.8, green squares: z = 0.8 — 1.6, blue 
diamonds: z = 1.6 — 2.8. An indication of the typical 68% error 
is shown offset from the real data points. The line where r par t 
= Tfull is marked on the figure for reference and does not rep- 
resent a particular linear fit. The solid green/dashed blue lines 
represent 68% error bounds on the best-fitting trendline for the 
z = 0.8 — 1.6/z = 1.6 — 2.8 sources respectively. It can be seen 
that all sources lie within these boundaries, suggesting the sub- 
samples are consistent with each other. 



function of hard (2-10 keV) X-ray luminosity which is plot- 
ted in Figure 1121 We use the hard band since this energy 
range is less affected by intrinsic absorption. 

We calculate Spearman's rank correlation coefficient, p, 
for both the entire population of 734 sources, and the subset 
of 650 sources whose best fit is a simple power law model 
(po). Considering the entire population we find a negative ~ 
8<7 correlation of T with log Lx (p = -0.24, prob = 4x 10~ n ). 
Considering only the sources fit by po, we find a tighter 
negative ~ 8a correlation (p — -0.29, prob = 2 x 10~ 14 ). 
These correlations imply that it is the sources with higher 
hard X-ray luminosities that tend to have natter F values. 

As previously discussed, the higher redshift objects in 
the sample may be present due to pre-selection. Since these 
sources are also likely to be the sources with the highest lu- 
minosities, we test whether their exclusion from the sample 
has any effect on the correlation coefficients. We find that 
when only considering sources with z < 2.6, the correlation 
of r with Lx remains present at 8a. 

Figure [13] shows a binned version of Figure 1121 with 
the top plot including sources fit with all models, and the 
bottom plot including just those fit with the simple power 
law model. Each bin size is approximately log Lx = 0.5, and 
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Figure 12. This plot includes the 734 sources with H0> 1%. The 
different colours/symbols correspond to the best fitting model; 
red circles are po, green squares are apo, blue diamonds are 
po+bb and black triangles are apo+bb. The error bars shown 
at the top of the plot correspond to the average of each of the 1<t 
errors on the individual sources. 



includes at least 18 sources. A weighted mean value of F is 
calculated for each bin using Equation [1] as before. 

As in the previous section, we calculate linear trends to 
the data using \ 2 minimization and also consider whether 
the data are consistent with a (r) that is constant with Lx ■ 
These trends along with the null hypothesis probabilities for 
each model can be found in Table [2] 

Considering sources best fit with all types of spectral 
model, the linear trend gives a gradient only consistent with 
zero at 3a. However, a constant F also gives an acceptable 
fit to the data meaning that a non-evolving V model cannot 
be ruled out. When only the sources fit with a simple power 
law model are considered, the non-evolving T model does 
not provide a good fit to the data and the linear trend gives 
a gradient inconsistent with zero. This suggests that for this 
sub-set of sources, a significant trend of T with Lx is present. 

We also investigate whether the trends we observe are 
still present in the data when we restrict the sample to 
sources above certain count levels. For the un-binned data, 
a significant correlation is still present between Lx and V 
even when all the sources with < 1000 counts have been re- 
moved, although the significance level does fall. The binned 
data shows similar results for all count levels i.e. a trend 
that is only consistent with zero at 3a (all models) and 5a 
(po models). 

We test whether the correlation we observe is still 
present when the RLQs are removed from the sample, since 
it has been suggested that they have flatter spectra than 
RQQ. A significant trend between F and Lx is still seen 
(5a). 

Since we find the stronger of the two correlations is 
between F and Lx, we suggest that it is luminosity that is 
the dominant variable, rather than redshift. We expect such 
a correlation from a consideration of the spectral shapes and 
how the luminosity is calculated. For example, a source with 
a steep spectral slope will have less emission over the hard 
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Figure 13. These plots show the dependence of weighted T val- 
ues on hard (2-10 keV) luminosity. The sources are binned in 
approximately log Lx = 0.5 bins, although the lower and upper 
bins do include sources slightly below 43.0 and above 46.0 re- 
spectively. The solid horizontal lines indicate the weighted mean 
values of each sample and the dotted lines show the linear fit to 
the data. The top plot includes data from 734 sources with good 
fits using models with freely varying T values and a linear fit of 
(r) = (-0.06±0.02)logL x + (4.6±l.l). The bottom plot includes 
data from 650 sources with good fits using the simple power law 
model and a linear fit of (r) = (-0.15 ± 0.03)logL x + (8.5 ± 1.2). 
The y error bars are standard errors on the mean and the x errors 
indicate the width of the bin. 



X-ray range, compared with a source with a flat spectral 
slope. 



3.4 Radio Properties 

The radio loudness values (Rl = Fr/Fo) of our 
sources were determ ined using radio fluxes from FIRST 
( Becker et al.l Il995l) and the optical fluxes from SDSS 
I Adelman-McCarthvl 120071 ') as pre viously described i n Sec - 
tion [2] Using the classification of iKellermann et all (|l989h 
where Rl > 10 indicates the source is radio loud, 552 sources 
were confirmed as radio quiet and 75 as radio loud. The 
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properties of these sub-samples are investigated in this sec- 
tion. 

The distributions of absolute i band magnitudes for the 
RQQ and RLQ sub-samples are found to not be significantly 
different (KS significance = 0.39). In addition the redshift 
distributions are not significantly different (KS significance 
= 0.87). The X-ray luminosity distributions for RQQ and 
RLQ presented in Figure [14] (top) are found to be signif- 
icantly different (KS significance = 8 x 10~ 9 ), with RLQ 
having higher 2- 10 keV X-ray luminosities. Since the red- 
shift distributions are the same, this is likely due to a phys- 
ical difference in the sources rather than an observational 
bias from the sample being flux limited. This suggests that 
an additional mechanism of producing X-rays is present, 
which increases the total X-ray emission. This ide a has been 
previously sug gested by iReeves fc Turned |2000j ). The ad- 
ditional component is perhaps related to the jet and may 
be Synchtrotron Self Compton (SSC) scattering where the 
radio photons produced as synchrotron emission from elec- 
trons spiralling in a magnetic field, undergo inverse Comp- 
ton scattering off the sa me electrons, up to X-ray energies 
jBand fc Grindlavlll986l ). 

It is thought that the synchrotron emission produced 
in RLQ will contaminate the power law produced by in- 
verse Compton scattering, making the X-ray spectrum ap- 
pear flatter. We therefore explore the differences in the V 
distribution of the sub-samples of RLQ and RQQ, as shown 
in Figure [14] (middle) and a KS test finds that they are 
significantly different (KS significance^ 3 x 10~ 5 ). When 
modelled with a single Gaussian, the RQQ distribution has 
a mean value of F = 2.01 ± 0.01, whereas the RLQ distri- 
bution has a lower mean value of F = 1.86 ± 0.02, however 
we note that the 2 populations cover a similar range of F 
values. Flatter Y values in RLQ are noted in the literature, 
although the actual values found here are steeper than those 
reported (e.g. I Reeves fc Turner|[2000l '). 

We also note that the dispersion of F values is much 
smaller for the RLQ than the RQQ. We find the intrinsic 
dispersions to be (T(y) = 0.31 ± 0.01 for the RQQ, which 
is similar to that for the entire original sample and a^ 
= 0.20 ± 0.02 for the RLQ, which is considerably lower. 
The typical 68% error on Y for the RLQ is 0.08, signifi- 
cantly lower than the intrinsic dispersion. Figure [14] (bot- 
tom) shows the best-fitting (r) and cr^ values and la, 2a 
and 3a confidence contours from the Gaussian fitting of the 
F distributions. This clearly shows that the Y distribution 
for RLQ has a lower dispersion than that for RQQ and the 
average (r) value is also lower. 

Figure [151 (top) plots the radio loudness parameter Rl 
against Y. The solid line shows the average Y value for the 
RLQ sample whilst the dashed line shows the best-fitting 
linear trend to the RLQ. We expect to see flatter Y val- 
ues for increasing levels of radio loudness as suggested from 
Figure Q3] (middle) and the idea that the power law will be- 
come contaminated by emission from the radio jet. We find 
the data is inconsistent with a constant F fit (p=0%), but an 
F-test shows that the linear fit is not statistically required. 
This is also the case when the combined RQ and RL popula- 
tions are considered. Since both the constant and linear fits 
give poor fits to the data due to the large spread in values, 
we also consider the data when it is binned by radio loud- 
ness. In this case the data is well fit by the constant Y model 
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Figure 14. These plots compare the Lx (top) and T (middle) 
distributions of the confirmed radio quiet (solid blue) and radio 
loud (dashed red) sources. The bottom figure shows the best- 
fitting (r) and a ij-) values (crosses) and la, 2<j and 3<r confidence 
contours from the Gaussian fitting of the T distributions for the 
RQQ and RLQ sub-samples. 
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(p=15.3%) and an F-test does not require the linear model. 
This is again also true for the combined RQ and RL pop- 
ulations and therefore we do not report a significant trend 
between Rl and T. 

Figure [15] (bottom) plots the radio loudness parameter 
Rl against Lx- The solid line shows the average Lx value 
for the RLQ and the dashed line shows the best-fitting lin- 
ear trend to the RLQ sample. The result from Figure [14] 
(top) suggests that the more radio loud objects also tend 
to be more X-ray luminous. The data is inconsistent with 
a constant luminosity model (p=0%), but a linear fit is not 
statistically required. We note that for the combined RQ 
and RL populations, the linear trend is required, but this 
model still does not provide a good fit to the data (p=0%). 
When we consider binned data, the constant log Lx model 
is still a poor fit to the RLQ (p = 5 x 10 -11 ), but the linear 
model is now statistically required. However, the gradient 
of the linear trend is 0.17 ± 0.08 which is only significant at 
~ 2a. For the combined RQ and RL populations, the linear 
trend is not required for the binned data. 

3.5 Intrinsic Cold Absorption 

Whilst the majority of our sources were well fit by a sim- 
ple power law model, some sources were better fit with an 
absorbed power law which includes an additional intrinsic 
neutral absorption component at the redshift of the source 
(modelled by zphabs in XSPEC). The requirement of this 
component was tested by the F-test at 99% significance. 

In the sample of 734 sources which have H0>1%, 32 
sources required an additional absorption component. Set- 
ting the F-test confidence at the 99% level, requires that 
we take into account the presence of 1% of detections which 
could be spurious. Accounting for these gives 3.4% absorbed 
sources. This represents only a lower limit to the number of 
absorbed sources since the detection of absorption at the 
required significance is more difficult in spectra with lower 
numbers of counts. This value is lower than those in the 
literature, with many studies reporting absor bed fractions 
of ~ 10 %, albeit for low e r F-test significances dCorral et al 



201 ll: l Young et alj|2009l ; iMainieri et al.ll2007l ; iMateos et al 
2005bl lalL iMateos et al] (|201Ch uses the more stringent sig 



nificance threshold of 99% as we do here, but still find a 
greater absorbed fraction of ~ 8%. 

Both the literature results and ours found here are at 
odds with the s tandard orientation based Unified Model 
|Antonuccilll993l ). which does not predict sources that are 
optically classified as type 1 to require an absorbed power 
law model in X-rays. The typical levels of absorption in our 
sources range from 10 21 cm -2 to 10 23 cm -2 , far higher than is 
expected for an AGN defined as type 1. Figure [161 shows the 
distribution of N} I ltr values for all the absorbed sources in 
our sample. A visual inspection of the optical spectra of the 
absorbed sources hints that some may be of an intermediate 
type (e.g. Sy 1.5), and others may be BALQSOs, offering 
a possible explanation for their apparent high levels of ab- 
sorption. An analysis of the optical reddening properties of 
the sample will be discussed in Hutton et al. (in prep). 

In Figure [17] (top) we investigate the dependence of the 
intrinsic column densities of the absorbed sources, measured 
in the rest frame of the source, on redshift. Although we see 
a slight trend for higher absorption levels in higher redshift 
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Figure 15. This plots the radio loudness, Rl, as a function of the 
best-fitting V values (top) and log Lx (bottom). The RLQ arc 
shown as red squares and the errors on T correspond to the 68% 
confidence level. The RQQ are shown as blue circles, for which 
the errors on T have been ommited for reasons of clarity. The 
solid lines show the average T (1.80) and log Lx (44.88) values 
for the RLQ sample, which represent the non-evolving models. 
The dashed lines show the best-fitting linear trends to the RLQ: 
T = (-0.022±0.006)logR c + (1.85±0.02) and \og(L x ) = (0.17 ± 
0.002)logi?. L + (44.5 ± 0.003). 



sources, this can be explained as an observational bias, since 
sources with low amounts of absorption will not be detected 
at higher redshifts due to their absorption signature falling 
outside of the EPIC bandpass. We do not detect absorption 
in any sources above z = 3 where the lower energy limit of 
the XMM-Newton bandpass corresponds to a source energy 
of ~ 2 keV. At z = 2, where this limit falls to ~ 1.5 keV, 
some absorbed sources become detectable, but only at levels 
high enough to create an absorption signature that falls into 
the observed bandpass. Considering correlation coefficients, 
we find no significant trend of intrinsic column density with 
redshift (prob= 0.02). 

We also consider the relationship between N l ff tr and 
hard X-ray luminosity in Figure [17] (bottom). We consider 
the X-ray luminosity in the hard band (2-10 keV) rather 
than the soft band (0.5-2.0 keV) as this is the energy range 
at which any intrinsic absorption present in the source would 
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Figure 16. The distribution of intrinsic absorbing column den- 
sities for the 32 sources in the sample that require the additional 
absorption component. 



suppress the emission. There is no correlation between jV™ 
and Lx which is confirmed by correlation coefficients (prob= 
1.0). 

We investigate how the fraction of absorbed sources 
changes with increasing luminosity and redshift. This can 
be seen in Figure [TH] and the fraction is found to be con- 
stant across the di fferent luminos ity and redshift bins. Pre- 
vious studies (e.g. iHasingerl 120081 ') which include both 'un- 
absorbed' type 1 and 'absorbed' type 2 objects find the 
percentage of absorbed sources decreases with increasing 
luminosity , whic h supports the Receding Torus model of 
iLawrenca (|l99ll ). This notes that the dust sublimation ra- 
dius increases as a function of luminosity and therefore the 
inner edge of the torus will be further away for higher lumi- 
nosity sources, resulting in less obscuration. Although our 
sample includes only type 1 objects, a small percentage of 
these have been shown to include intrinsic absorption, which 
might be expected to behave in a similar way to the ab- 
sorption found in type 2 objects. However, since we show 
the percentage of absorbed objects to remain constant with 
increasing luminosity, the receding torus model is not sup- 
ported and the absorption present in our type 1 sources must 
be of a different nature. 



3.6 Soft X-ray Excesses 

The sources were also fit to include a soft excess component, 
modelled with zbbody in XSPEC (see Section [OJ . We note 
that using a blackbody to model the soft excess component 
is not strictly physical, but provides a good representation of 
the spectral signatures and allows a simple comparison with 
other studies. We only search for this component in sources 
with X-ray spectra containing at least 100 counts. Using the 
F-test at 99% significance, we detect this component in 60 
of 680 sources. Accounting for spurious detections, we find 
evidence for a soft excess component in 7.8% of the sample. 
This fraction agrees with p revious studies of type 1 AGN 
(|Mateos et all l20ld . l2005aT ). We note again that this is a 
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Figure 17. The relations between the intrinsic absorbing column 
density and redshift (top) and 2-10 keV X-ray luminosity (bot- 
tom). Green squares indicate the sources best fit with the apo 
model and black triangles indicate the apo+bb model. The error 
bars correspond to the 68% confidence range. 



lower limit of the percentage of sources we expect to have a 
soft excess component; this value could be as high at 80%. 

In the automatic fitting process, XSPEC is left to fit the 
po+bb model, which has 5 free parameters, with few con- 
straints. In 7 cases a flat power law and a large blackbody 
component gives a lower value of \ 2 , despite this fit being 
less physically realistic. These sources were re-fit over the 
0.5-12.0 keV detector frame energy range with the po+bb 
model, keeping the V value fixed to that obtained from a fit 
to the data over the 2-10 keV rest frame energy range using 
just a simple power law. V and kT values more consistent 
with the other sources are then obtained. An F-test between 
the manual fits and the simple power law fits is not appro- 
priate since the new fit utilizes a result already obtained 
from the data. It is therefore assumed that the blackbody 
component is still statistically required in the fit, based on 
the original classification. 

The mean temperature of the soft excess component 
when fitted with a blackbody is (kT) = 0.17±0.09 keV, ob- 
tained from modelling the distribution with a single Gaus- 



16 Scott et al. 



" 20 



w 15 



• 



42 



43 



44 45 46 

loq U 



• 5 



' 



L. 
0.0 



0.5 1.0 1.5 2.0 2.5 3.0 

Redshift 



Figure 18. These figures show how the percentage of sources that 
require an intrinsic absorption component in their best-fitting 
model varies for different luminosity and redshift bins. The x er- 
ror bars indicate the width of the bins, and the y error bars come 
from Poissonian counting statistics. 
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Figure 19. The distribution of kT values for the 56 sources in 
the sample that require the additional soft excess component and 
whose best fit provides an acceptable fit to the data. 



sian. This is slight l y higher than results in the literature 



llMateos et alJ20ld:lwinter et al.ll2009l : iMateos et ai1l2005bl ; 
iGierlinski fc DondbOO^ although consistent within errors. 
The range in temperatures is observed to be rather small, 
with 68% of the sources having kT values in the range 0.1- 
0.25 keV. The full distribution of kT values can be seen in 
Figure [H 

We consider possible correlations between the soft ex- 
cess temperature log(kT), and redshift, hard X-ray luminos- 
ity, the luminosity of the blackbody component (Lbb) and 
the black hole mass. The correlation values and best-fitting 
linear trends can be found in Table [3] We find no signifi- 
cant correla tion with z, which is consistent with previous 
results (e.g. IMateos et al.|[201ol '). Figure [20] (top left) shows 



the presence of some soft excess components detected even 
in sources at z > 1. 

The standard explanation for the soft excess, and the 
motivation for modelling it with a blackbody is that it is 
thermal emission from the inner accretion disc. The stan- 
dard description predicts characteristic temperatures of the 
order kT ~ 0.01 keV for a 1O 7 M0 black hole, as determined 



from T, 



disc 



M 



-1/4 



(L/Li 



a/4 



(jShakura fc Sunvaevl 



1 19731 ). We note that our average kT value is too high to be 
consistent with the standard disc origin. This scenario also 
suggests that the temperatures should show variations with 
luminosity and black hole mass, with positive and negative 
linear trends, respectively. 

We find a significant correlation with Lx (see Figure [2H1 
bottom left), but no significant correlation with Lbb- A 3a 
correlation is found between kT and Msfl[f], but the best- 
fitting linear trend shows a positive, rather than negative 
gradie nt as would be expec ted. This result disagrees with 
that of iBianchi et alj (|2009h who report the temperature of 
the blackbody component used to model the soft excess to 
be unrelated to the mass of the black hole. 

Previou s studies have reported no trend of kT with 
luminosity (|Bianchi et al.1 12009| ; Iwinter et al.1 12009| ) This 
pro blem was first highlig hted by Gierlihski fc Done] (|2004T ) 
and lCrummv et all (|2006t l. which led to the development of 
new theories, moving away from the accretion disc origin. 
They invoke the idea that since the soft excesses display 
a narrow range of kT values, the physical process behind 
their origin is likely to be one of constant energy, such as 
atomic transitions. At ~ 0.7 keV, lines and edges of ion- 
ized O VII and O VIII produce a strong jump in opacity. 
If these features were smeared by high velocities or gravita- 
tional redshifts, as found close to a black hole, a featureless 
soft excess would be produced. Therefore the soft excess 
component could be an artefa ct of ionized absorption, likely 
from an outflowing disc wind (|Gierliriski fc Donell2004f ). Al- 
ternatively, the partially ionized material could be out of the 
line of sight and seen via reflection. In this scenario, the hard 
X-ray source needs to be at a small height above the disc so 
that light bending can illuminate the inner disc. In this case 
the soft excess component is an art efact of relativistical ly 
blurred photoionized disc reflection l|Ross fc Fabianll2005l ). 

In Figure [2T1 we plot black hole mass against the black- 
body luminosity and find a ~ 3cr correlation. This suggests 
that regardless of the origin of the soft excess, its luminosity 
is likely related to the mass of the black hole. 

Figure [22] shows a strong (~ 8a) correlation between 
the luminosity of the blackbody component and the lumi- 
nosity of the power law component with a best-fitting trend 
of logLbb = (0.78±0.02)logi PL + (9.2±1.0). This indicates 
a link between the energy production in these 2 components 
and supports the idea of inverse Compton scattering where 
the low energy photons from the disk are up-scattered to 
X-ray energies creating the power law emission. However it 
could also imply a reverse process in which the power law 
emission is re-processed into the lower energy soft excess. 
This correlation has been previously reported in 



3 For 29 of the sources that required a blackbody component in 
the spectra l fit, w e also have a black hole mass estimate from the 
Sh en et al.1 1 |2008|) catalogue (See Section [32] for details) 
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Table 3. The Spearman rank correlation coefficients and significances for correlations between the characteristic blackbody temper- 
ature used to model the soft excess, kT, and the parameters z, Lx, £bi> and Mbh- The best-fitting linear trends were determined 
by x minimization. 
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Figure 20. The relations between the characteristic soft excess temperature, kT, and redshift (top left), hard X-ray luminosity (bottom 
left), blackbody luminosity (top right) and black hole mass (bottom left). The error bars correspond to the 68% confidence range. Only 
sources for which H0>1% are included on these plots. Sources indicated by the filled black circles or upper limits are the ones fit with a 
manual po+bb model as mentioned in the text. These sources are not included in the determination of any correlations or linear trends, 
however they are shown to be consistent with the rest of the sample. 



IWinter et all (|2009h . where they find L PL = (0.79 ± 
0.14)Li,i) + (9 . 34 dz 6.04). A similar correlation is seen in 
iMateos et all (|2010l) where the power law luminosity over 
2-10 keV against the blackbody luminosity over 0.5-2 keV 
is plotted. The best-fitting linear function is found to be 
logLti, = (1.2±0.1)logLp L + (9.2±5.3). It is suggested that 
the lack of sources in the bottom right of the plot is a selec- 
tion effect due to the fact that a low luminosity blackbody 
component will be difficult to detect over a high luminos- 



ity power law component and the lack of sources in the top 
left of the plot is due to an apparent upper limit to the soft 
excesses allowed in AGN. For comparison, the linear trend 
best fit to our data when the errors are not considered is 
logLw, = (0.80 ± 0.09)logL PL + (8.3 ±4.1). These linear re- 
lationships each with gradients close to 1 imply that the 
fraction of converted emission is constant for sources of all 
luminosities. However all 3 data sets show that in the lumi- 
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Figure 21. This figure includes the 25 sources best fit with the 
po+bb model, for which we also have a black hole mass estimate 
and a determination of Lbb- The different colours correspond to 
the different emission line used in the mass determination (red - 
H/3, green - Mgll). It shows the dependence of the black hole mass 
on the luminosity of the blackbody component, determined over 
the 0.5-10.0 keV range. The errors shown on L(,j, correspond to 
the 68% confidence range. The black circles indicate the sources 
fit with the manual po+bb model. 



nosity range 10 — 10 erg s , the luminosity in the power 
law is greater than that in the blackbody component. 



3.7 Accretion Properties 

Since the X-ray emission is believed to originate from the 
accretion process, relationships between the X-ray spectral 
properties and the physical properties of the central engine 
are expected. In this section we consider correlations be- 
tween both the broad-band continuum and the soft excess 
properties with Eddington ratio. We also investigate any dif- 
ferences in Eddington ratio between the radio loud and radio 
quiet sub-samples. 

Black hole masses can be determined for type 1 AGN by 
assuming that the broad line region is gravitationally bound 
to the central black hole potential and that the broad line 
clouds are virialized. By measuring the width of broad emis- 
sion lines produced by these clouds their velocity (V 9 ) can be 
estimated. The radius of the broad line region (Rblr) can 
be estimated from a linear correlation with the monochro- 
matic optical luminosity. The mass of the black hole can then 
be determined fro m the equation, Mbh = RblrVqG -1 
<|Kaspi et al.ll2000h . 

Here we use the IShen et al.l (|2008l ) catalogue of black 
hole masses estimated using this technique, which includes 
~ 58, 600 quasars from the SDSS DR5 quasar catalogue. 
The mass is estimated using a different broad line, depend- 
ing upon the redshift of the source. H/3 is used for z < 0.7, 
Mgll for 0.7 < z < 1.9 and CIV for z > 1.9. The catalogue 
also includes bolometric luminosities estimated from optical 
monochromatic luminosities (at 5100A for H/3, 3000A for 
Mgll and 1350A for CIV), derive d usin g bolometric correc- 
tion factors from iRichards et al. | (|2006h . Of the 734 sources 
(with HO > 1%), 79% have a black hole mass estimate avail- 



ed 
cn 
CD 



O 

c 

E 

_J 

>N 

X3 
O 
X3 

u 
m 



10 4 



10 4 



10 4 



10* 



L„ 


= 100% L PL 


U, 


= 50% L, L 


L„ 


= 10% L PL 




10 4 

1 42 1 Q 43 1Q 44 1Q 45 ^46 

Power Law Luminosity (erg s _1 ) 

Figure 22. This plot compares the luminosity included in each 
of the separate components of the po+bb fits, calculated over 
the 0.5-10.0 keV range. The linear trendline is logLbb = (0.78 ± 
0.02)logLp£ + (9.2 ± 1.0). The errors shown in these plots cor- 
respond to the 68% confidence range. The black circles indicate 
the sources fit with the manual po+bb model and are consistent 
with the rest of the sources. The lines indicate the locations of 
Lbb = (solid), L^b = 50%Lpi (dashed) and Lbb = 10%Lp£ 

(dotted) but do not represent any particular fit to the data. 



able in the lShen et al.l (|2008l ) catalogue. This gives us a sam- 
ple of 581 sources. 

The Eddington luminosity is calculated using Equa- 
tion [3] The bolometric luminosity is calculated from the 
observed X-ray luminosity in the 2-10 keV band, Lx, 
using the lum i nosity dependent bolometric correction of 
iMarconi et al.l |2004h described by Equation [4] where C = 
logLtoi — 12 and Lboi is in units of L©. The Eddington ratio 
is then defined as Lboij L^dd- 



AnGMm p c 

L<Edd = 



1.3 x 10 3S — ergs" 1 

M© 



lag(L bo i/L x ) = 1.54 + 0.24£ + 0.012£ 2 - 0.0015£ 3 



(3) 
(4) 



In Figure[23]we plot the relationship between T and the 
Eddington ratio, separated depending on which broad line 
was used to estimate the black hole mass. We find a positive 
correlation for the H/3 data (Spearman Rank correlation co- 
efficient, p — 0.26, 2.8a), no correlation for the Mgll data 
and a negative correlation for the CIV data (p = —0.32, 
2.9<t), however the reliability of the CIV line as a mass es - 
timator has been widely questioned (e .g. IShen et alJ|2008h . 
This differs slightly with the results of lRisaliti et al.l ( 2009f ) 
who find a highly significant, positive correlation with esti- 
mates based on the H/3 line, a significant, although weaker 
positive correlation with Mgll and no significant correla- 
tion with CIV. We note that in their analy s is the y use the 
bolometric luminosity given by IShen et al.1 (|2008l ) and use 
the power law slope determined in the 2-10 keV rest frame 
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band which may explain these differences. Previous studies 
of PG Quasars with XMM-Newton have repo rted a positive 
correlation between F and Ed dington rate (IPorcmet et al.l 
12004 IPiconcelli et al.ll2005l ) and lShemmer et al.l (|2006t ) show 
that this is the dominant correlation, as opposed to the 
known relation between accretion rate and line width. We 
note that whilst we use a luminosity dependent bolometric 
correction to obtain our bolometric luminosity estimates, 
recent works have indicated that the bolometric correc- 
tion factor may be depende nt on the Eddington ratio itself 
jVasudevan fc Fabia J2007h . 

The previous section indicated a strong relationship be- 
tween the luminosities in the blackbody component and the 
power law component. Since the properties of the power law 
component may correlate with the Eddington ratio, we also 
investigate the dependence of the blackbody properties on 
Eddington ratio. We consider the dependence between kT, 
blackbody luminosity and black hole mass on the Edding- 
ton ratio. No significant correlations are found (kT = 0.8a, 
L bb = 2.4<j, M BH = 0.3<t). 

In order to investigate a possible link between radio 
jet production and accretion rate, we have compared the 
distribution of Eddington ratios for the confirmed RLQ and 
RQQ sub-samples and investigated a possible trend between 
radio loudness and Eddington ratio. We consider both an X- 
ray derived Eddington ratio, (as previously described), and 
an optically derived Ed dington ratio usin g the bolometric 
luminosities derived by IShen et alj (|2008l ) (also described 
above). The distributions of Eddington ratios are plotted in 
Figure [24] (top) and the KS significances are listed in Ta- 
ble!?] Figure [5?] (bottom) shows the correlation between ra- 
dio loudness, Rl and Eddington ratio, with the best-fitting 
linear trends listed in Table [4] 

We find that the distribution of X-ray derived Edding- 
ton ratios is significantly different between the RQQ and 
RLQ sub-samples, with RLQ having higher Eddington ra- 
tios. However, when the optically derived Eddington ratios 
are compared, no significant difference is found. Similarly, 
a significant trend towards higher X-ray derived Eddington 
ratios in sources with a higher radio loudness is suggested. 
This trend is not found when optically derived Eddington 
ratios are used. 

We suggest that the differences found between using X- 
ray derived and optically derived Eddington ratios may be 
an artefact of the additional X-ray component suggested to 
be present in radio loud quasars (see Section [3. 4p . Whereas 
the optical emission of RQQs and RLQs is similar (the dis- 
tributions of Mi were found to not be significantly different), 
the X-ray luminosity of the radio loud sources was found to 
be higher, which results in relatively higher X-ray derived 
Eddington ratios for these sources. 
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Figure 23. This figure shows the relationship between T and 
the Eddington ratio, where the latter is based on mass estimates 
obtained using different broad emission lines. Top: 117 sources 
for which the H/3 line was used. Middle: 385 sources for which 
the Mgll line was used. Bottom: 79 sources for which the CIV 
line was used. The error bars correspond to 68% errors. Sources 
indicated with the black dots are RLQ. 



4 DISCUSSION 

Much of the information gained to date on the X-ray emis- 
sion mechanisms of AGN and their interplay with the sur- 
rounding environment has come from detailed studies of 
brig ht AGN, which reve al significant spectral complexities 
(e.g. iPounds et aT]|2004h . The nature of our sample here 
means that despite dealing with lower count spectra which 
allow only relatively simple models to be fit, we have a suf- 



ficiently large sample with which to investigate source-to- 
source variations and the range in spectral properties. 

At the simplest level the sample is well described by a 
simple power law model. There is significant intrinsic dis- 
persion in the power law indices describing the broad-band 
X-ray continuum, with a small fraction of objects which are 
either extremely flat or extremely steep. We emphasize that 
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Table 4. Listed in the first column arc the KS significance levels found when comparing the distribution of Eddington ratios 
(derived from both X-ray and optical bolometric luminosities) for RQ and RL sub-samples. Also listed are the best-fitting linear 
trends between radio loudness, R^, and Eddington Ratio, for the combined RQ and RL populations, and just the RLQ on their 
own. 
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Figure 24. The top figures show the Eddington ratio distributions of the confirmed radio quiet (solid blue) and radio loud (dashed red) 
sources. The bottom figures show the trend of R^ with the Eddington ratio with the best-fitting trend to the RLQ shown by the dashed 
line. The figures on the left use Eddington ra tios defined from the observed 2-10 keV X-ray luminosity and the luminosity dependent 
bolometric correction of llVfajcoiuet aL I J2004) . The figures on the right use Edding ton ratios defined from monochromatic luminosities 
in the optical and UV b-y lSrien*^t**aT r^oosD , and bolometric correction factors from lRichards et al. (2006). 



for a number of these extreme objects we have sufficient 
counts to rule out effects caused by the non-inclusion of ad- 
ditional spectral complexity. The distribution is reasonably 
symmetric about the mean, suggesting that the underlying 
electron acceleration mechanism and/or Comptonisation are 



subject to a rand om scatter, as opposed to a limiting cas- 
cade mechanism (ISvensso ni l 19941 ) which might be expected 
to produce a larger asymetry in the distribution of values. 

Despite evid ence that nearby AGN display rapidly vari- 
able spectra (e.g. IVaughan et al.|[201ll ) , we assume constant 
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X-ray spectra for our so urces. However, a Prin cipal Com- 
ponent Analysis (PCA; iFrancis fc Wills! Il999l ) on a long 
XMM-Newton observation of MCG-6-30-15, found that 96% 
of the variability is due to the continuum component, the 
spectrum of which is well fit with a simple power law of 
photon index F = 2.20 ± 0.05. The spectral variations are 
thought to be caused by the normalization of this power 
law component changing relative to a constant component, 
rather than changes in the underlying power law index 
jVaughan fc Fabianlliool ). The value of this is found to be 
similar to the average F value of our sample. 

The concept of Black Hole Unification says that since 
AGN and Galactic Black Hole Binaries (GBHBs) are both 
examples of accreting black holes, they should be expected 
to display similar characteristics, in which timescales sim- 
ply scale with mass. In particular, GBHBs exist in differ- 
ent spectral states corresponding to different accretion rates 
with transitions between states occuring on the order of 
weeks. For AGN this timescale is ~ 10 years, meaning that 
such transitions are not observed for individual AGN. An 
analogy is made between the spectral states of GBHBs and 
the different types of AGN with some RLQ perhaps being 
objects transitionin g to the high/soft state via the jet line 
jFender et al.ll2006l ). We expect the majority of objects in 
our sample to be in the high/soft state due to selection bi- 
ases. 

The most significant factor to correlate with spectral 
index is the intrinsic X-ray luminosity of the object. This 
finding is similar to that of I Young et alj (|2009h who report 
no correlation with redshift, but a strong correlation with 
luminosity. They report strong correlations of V with rest- 
frame monochromatic luminosities at 0.7, 2.0, 10.0 and 20.0 
keV, with the strength of the correlation increasing towards 
higher energies. However, since the monochromatic luminos- 
ity determination depends directly on the F value, a correla- 
tion between the two quantities is expected. As a flat source 
has more emission at higher energies than a steep source, 
it is the natter sources that will have the higher monochro- 
matic luminosities at these energies. This effect will decrease 
between hokev an d hkev, until a pivot point at which the 
steeper sources start to have more emission at the energy at 
which the monochromatic luminosity is calculated, than the 
flat sources. 

Calculating the Eddington ratio requires an estimate 
of the bolometric luminosity of the source. This can be in- 
ferred from either the observed X-ray or optical luminosity 
and we obtain different estimates depending on which wave- 
band is used. In particular we show that RQQ and RLQ 
have similar optically derived Eddington ratios but higher 
X-ray derived Eddington ratios, likely due to the generally 
higher 2-10 keV X-ray luminosities seen in RLQ, on which 
the ratio is based. Whilst we note that an inherent bias due 
to the RLQ is present, we see a trend for an increased Ed- 
dington ratio with increasing radio loudness, which might 
argue for models in which jet formation and hence radio 
emission depend on the accretion mode/Eddington ratio in 
analogy with galactic stellar mass black holes. We also see 
some evidence for the correlation of spectral index with ra- 
dio loudness. In particular we show that the RLQs have on 
average flatter power law slopes which may be related to 
the synchrotron emission from a jet whose importance de- 
pends on the radio luminosity. Therefore for objects with 



higher Eddington ratios we expect a higher radio loudness 
and from this a flatter power law spectrum. There is how- 
ever conflicting evidence for a correlation of the power law 
index with the Eddington ratio. 



In agreement wit h previous works ( Porquet et al.ll2004l; 



Piconcelli et al. | 120051 ; [Sh emmer et al.l 



2006; iRisaliti et ail 



20091) we find a, positive correlation between F and the Ed- 



dington ratio. However, this is only true when the black 
hole mass estimates are determined using the width of the 
broad H/3 line; if the Mgll line is used, no correlation is seen, 
and the correlation appears to switch direction when using 
the CIV line. The virial method for determining black hole 
masses is calibrated using results from reverberation map- 
ping studies. These are only widely available for low redshift 
objects, giving a well calibrated method using the H/3 broad 
line but a less reliable calibration for the Mgll and CIV lines 
whic h are simply scale d to match the results from the H/3 
line. IShen et al.l (|2008l ), from whom the black hole mass es- 
timates are taken, note that the estimates determined from 
the H/3 and Mgll lines tend to agree with each other, but 
the masses estimated from the CIV line do not. The CIV 
lines can be severely affected by a disc wind component and 
therefore may be asymetric and blue shifted, leading to an 
overestimate of the width of the line and a FWHM which 
is not representative of the virial velocity of the broad line 
region. However, some authors have reported consistency be- 
tween blac k hole masses determined with the H /3, Mgll and 
CIV lines (|Assef et al.ll20ld ; lGreene et al.ll2010T ). Whilst the 
different correlations we see may be due to a calibration ef- 
fect, we cannot rule out the possibility of a hidden physical 
effect. By necessity different broad lines are used for sources 
at different redshifts and the CIV line is used for sources at 
high redshifts and likely high X-ray luminosities. Therefore 
the difference in these parameters may be the true driver 
behind the different trends we see with Eddington ratio. 

Over recent years there has been much debate about the 
origin of the soft excess detected below ~ 0.5 keV in AGN. 
Initially per ceived as simply blackbo dy emission from the ac- 
cretion disc (|Turner fc Pound3ll989l ') the effective blackbody 
temperatures measured are significantly higher than accre- 
tion disc models allow. Alternative models have been pro- 
posed, including Comp ton upscattering of the disc emission 
|Kawaguchi et al 1l200ll). reflection o f the power law emission 
from the disc ( Crummy et al.|[2006l ). an d absorption of the 
powe r law emission by ionised material (|Gierliriski fc Donel 
I2004D . Our results, based on fitting this component with a 
blackbody, argue against, but do not rule out, these absorp- 
tion models, for which a correlation between the apparent 
temperature of the blackbody and fractional luminosity of 
the apparent excess might be expected. We find a strong 
correlation between the luminosities of the soft excess and 
the power law components, with the soft excess luminosi- 
ties being typically 10-50% that of the power law. This is 
easier to explain with an emission process such as Compton 
upscattering. 



5 CONCLUSIONS 

In this paper the X-ray spectra of a sample of 761 type 1 
AGN from a cross-correlation of the SDSS DR5 quasar cat- 
alogue and 2XMMi have been analysed. The main conclu- 
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sions from this analysis of the largest sample of AGN X-ray 
spectra are summarized below: 

(i) When the distribution of best-fitting power law slopes 
is modelled with a single Gaussian as done in previous works, 
the best-fitting values are found to be (F) = 1.99±0.01 with 
an dispersion of <j(r> = 0.30 ± 0.01. The typical 68% error 
on a source is AF = 0.13, which is much lower than the 
intrinsic dispersion of values. The distribution is better fit 
with a combination of 2 Gaussians due to the presence of 
extreme sources with flat or steep F values. 

(ii) A marginal trend (~ 3<j) for flatter F values in higher 
redshift sources is found, but only when just those sources 
best fit with a simple power law model are used. This appar- 
ent flattening does not appear to be due to the presence of 
a reflection component in the spectrum at higher energies. 

(iii) A significant trend (~ 5a) for flatter F values in 
higher X-ray luminosity sources is found. 

(iv) The radio loud quasars are found to have higher X- 
ray luminosities compared to the radio quiet quasars. We 
suggest this may be due to an additional X-ray component, 
perhaps related to a jet, which also causes RLQ to show flat- 
ter r values and higher Eddington ratios, if Lbol is calculated 
from the X-ray luminosity. 

(v) We find evidence for intrinsic neutral absorption in 
3.4% of our sample. However, our simulations suggest that 
the true fraction of absorbed type 1 sources could be as high 
as ~ 25%. The levels of absorption (N H = 10 21 -10 23 cm" 2 ) 
are typically much higher than the optical extinction of the 
broad-lines and those expected for type 1 objects. No trend 
is found between N}? tr and redshift or X-ray luminosity. 
The fraction of absorbed sources does not appear to vary 
with luminosity or redshift. 

(vi) We find evidence for a soft excess component in ~ 8% 
of our sample, but simulations suggest that up to 80% of 
our sources may have a soft excess. The average blackbody 
temperature used to model the soft excess is (kT) = 0.17 ± 
0.09 keV. This temperature is shown to correlate with hard 
luminosity. 

(vii) A strong correlation (~ 8a) is found between the 
blackbody luminosity and power law luminosity for those 
sources best fit with the po+bb model with the soft excess 
luminosity being ~ 10 — 50% that of the power law. 

(viii) A positive correlation (p = 0.26, 2.8<r) of F with the 
X-ray derived Eddington ratio is found when the broad H/3 
line is used in the determination of the black hole mass. No 
such trend is found when the Mgll line is used and a negative 
correlation (p = —0.32, 2.9<j) is found when the CIV line is 
used. The blackbody properties kT, Ltb and Mbh are not 
found to correlate with Eddington ratio. 
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